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ABSTRACT 


Doppler  spectra  of  radar  sea-echo  were  investigated  theoretically 
and  experimentally.  Electromagnetic  scattering  models  were  devel¬ 
oped  with  composite  rough  surfaces,  i.e.,  Eragg  resonant  water  waves 
superimposed  cn  a  carrier  water  wave;  and  basic  concepts  in  hydro¬ 
dynamics  were  developed  for  a  gravity  wave  of  finite  height.  The  ideas 
are  carried  over  to  the  sea,  where  the  assumption  of  a  fully  developed 
sea  yields  wave  height  expressions  of  the  form  ~//°  4  1  0  0  6  for  the 
width  of  the  spectra  of  radar  sea-echo,  H  being  the  height  of  the  carrier 
water  wave.  The  simple  models  yield  the  order  of  magnitude  and  the 
wave  height  dependence  of  the  spectra  width,  but  they  fail  to  give  its 
observed  polarization  dependence. 

New  experimental  results  on  spectra  of  radar  sea-echo  from  the 
Four -Frequency  Radar  (428  MHz,  1228  MHz,  4455  MHz,  and  8910  MHz), 
for  sea  conditions  which  ranged  from  1  to  2  ft  waves  and  1  to  2  knot 
winds  to  26-ft  waves  and  46  to  48  knot  winds,  show  that  the  bandwidth 
of  radar  sea-echo  is  polarization  and  frequency  dependent,  and  the  dif¬ 
ferential  doppler,  between  horizontal  and  vertical  polarization,  is  fre¬ 
quency  and  depression  angle  dependent.  Both  the  "noisiness''  of  the 
spectra  and  the  greater  spectral  width  of  the  horizontal  return  indicate 
spray  over  the  see  might  be  the  mechanism  responsible  for  the  depend¬ 
encies  observed. 
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STUDY  OF  DOPPLER  SPECTRA  OF  RADAR  SEA-ECHO 


INTRODUCTION 

During  the  past  years  the  scattering  mechanism  responsible  for  radar  sea-clutter 
has  been  investigated  in  great  detail.  As  a  consequence,  with  the  advent  of  recent  devel¬ 
opments  in  slightly  rough  and  composite  surface  scattering  theories  [1-5],  the  problem 
of  modeling  the  properties  of  average  radar  backscattering  cross  sections  seems  to  be 
well  under  control. 

To  further  improve  the  sea-clutter  model  it  is  desirable  to  have  additional  informa¬ 
tion  on  doppler  spectra  of  radar  sea-echo  to  gain  new  perspectives  on  the  electromagnetic 
scattering  processes  from  a  moving  air-sea  boundary.  A  review  of  previous  work  on 
doppler  shows  that  Crombie  l6j,  in  1955  at  13.56  MHz,  observed  that  the  mean  doppler 
frequency  of  radar  sea-echo  corresponded  to  the  phase  velocity  of  water  waves  which 
satisfied  the  Bragg  resonant  condition  [7],  i.e.,  water  waves  of  wavelength  half  the  radar 
wavelength,  for  grazing  incidence.  In  1958  Ament  et  al.  [8],  with  a  horizontaUy  polarized 
radar  at  22C  MHz,  noted  that  the  width  of  doppler  spectra  oi  sea-echo  was  comparable  to 
the  width  of  the  distribution  of  velocity  obtained  from  the  complete  energy  spectral  den¬ 
sity  of  the  water  waves  on  the  surface,  indicating  that  in  addition  to  the  resonant  scatter- 
ers,  other  water  waves  make  a  significant  contribution.  Hicks  et  al.  [9],  from  the  Uni¬ 
versity  of  Illinois,  using  horizontally  polarized  X-band,  concluded  that  radar  return  was 
produced  by  two  types  of  scatterers,  one  set  made  by  "unbroken"  water,  of  dimension 
one-half  the  wavelength  of  the  radar  frequency,  moving  at  the  orbital  velocity  (fluid  ve¬ 
locity)  of  the  gravity  waves  of  maximum  amplitude,  and  a  second  set  made  by  wind-blown 
spray.  The  half-power  width  of  the  spectra,  as  a  function  of  sea  conditions,  followed  the 
empirical  expression  11  Hs  ?m>  where  tf,  is  the  average  of  the  third  highest  waves  (sig¬ 
nificant  wa,Te  height)  and  Tr  is  the  period  of  the  gravity  wave  of  maximum  amplitude. 
Recently  Pidgeon  [10  ,  with  dual  polarized  C-band  coherent  radar,  for  grazing  angles  of 
0.3  to  8  degrees,  found  that  the  frequency  at  which  the  spectra  had  the  maximum  ampli¬ 
tude  was  polarization  dependent,  with  the  doppler  spectra  for  horizontal  polarization 
having  the  maximum  amplitude  at  a  higher  frequency  than  the  vertical  polarization. 

The  observations  at  220  MHz  and  higher  indicate  a  composite  type  of  surface  produc¬ 
ing  the  scattering,  i.e.,  resonant  water  waves  superimposed  on  gravity  waves  or  swells. 
Consequently,  in  the  sea,  the  transition  from  a  purely  slightly  rough  surface  to  a  compos¬ 
ite  rough  surface  (electromagnetkalJy)  must  occur  somewhere  between  13.56  MHz  and 
220  MHz. 

The  composite  surface  model  has  already  been  used  successfully  to  explain  polariza¬ 
tion  effects  ir.  radar  backscattering  cross-section  experiments  .4.5;,  and  its  applicability 
to  doppler  effects  will  now  oe  explored. 

Using  basic  concepts  in  hydrodynamics  to  develop  composite  surface  scattering 
models  for  gravity  waves  of  finite  height,  the  order  of  magnitude  of  orbital  motion  ef- 
fects  that  might  be  expected  in  the  spectia  of  radar  sea-erho  will  oe  estimated.  Finally, 
new  experimental  results  will  be  presented  on  doppler  spectra  obtained  from  the  phase 
channel  of  the  four-frequency  radar  of  the  Naval  Research  Laboratory,  both  to  verify  the 
models  and  show  that  the  sea-clutter  process  is  much  more  complex  than  the  earlier 
experiments  at  single  frequencies  have  suggested.  The  simple  models  account  only  for 
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orbital  iiiotion  effects  in  the  spectra,  while  other  mechanisms  seem  to  be  significant  in 
the  spectra  for  horizontal  polarization. 


SCATTERING  MODELS 

Some  basic  concepts  in  classical  hydrodynamics  will  be  introduced  first  in  order  to 
formulate  the  scattering  models.  In  the  sense  of  the  Stokes  expansion,  the  dynamics  of  a 
water  wave  of  finite  height  can  be  derived  from  the  dynamics  Of  a  sinusoidal  water  pro¬ 
file  [  11 }  (Fig.  1) 


£ ( x, t )  =  a  cos  A  (1) 

to  which  corresponds  a  velocity  potential  (infinite  depth) 

t)  =  ojak~le':z  sin  X  ,  (2) 

where  <u  is  the  angular  velocity  of  the  water  wave,  k  the  wave  number  of  the  water  wave, 
and  v  =  k  •  x  -  cot  the  phase  factor. 


To  the  first  order  of  the  wave  amplitude  a,  the  velocity  components,  u< 1  >  and  »< 15 
(Lagrangian  description)  are  identical  to  the  Eulerian  velocity  components,  dx/dt  and 

dx/dt: 


dx/dt  -  &4>/dx  -  uf  1  *  =  (jack*  cos  V  .  (3®) 

dz/dt  -  d<t>/dz  -  sin  X  .  (3b) 


Carrying  Phillips  procedure  to  the  second  order,  the  fluid  velocity  can  be  shown  to  be 

ut  =  u{  ♦  u^15  =  <uae*f  cos  X  ♦  ua3kelkl  .  (4a) 

wt  =  wj  1  >  ♦  wf  *>  =  <uaekt  sin  X  .  (4b) 

To  the  first  order  the  fluid  describes  circles,  and  to  the  second  order  a  drift  velocity 
appears  in  the  horizontal  component.  Kinsman  (12]  has  shown  that  in  the  Stokes  procedure 
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of  the  surface  profile,  *  -  £  is  no  longer  sinusoidal  but  to  the  third  order  tends  toward  a 
trochoid.  Thus,  a  trochoical  wave  profile  could  be  used  together  with  Eq.  (4)  for  the  fluid 
velocity  components  to  model  a  gravity  wave  in  the  sea.  Such  a  model  must  be  approxi¬ 
mate,  since  the  Stokes  expansion  has  been  shown  to  have  a  zero  radius  of  convergence 
and  the  Gerstner  wave,  which  has  a  pure  trochoidal  shape,  assumes  zero  drift,  which  is 
not  true  in  the  sea.  Kinsman  [  13]  notes  that  in  the  sea  the  waves  have  properties  of  both 
Stokes  and  Gerstner  wave.*  The  local  slope  for  a  trochoidal  profile  is  given  by 


r)£/dx  -  -  ak  sin  X/(  1  -  ak  cos  V)  ,  (5) 

which  for  waves  of  infinitesimal  amplitude  is  the  same  as  for  a  sinusoidal  profile. 

A  deterministic  type  of  composite  scattering  surface  can  be  constructed  using  a 
long-crested  carrier  water  wave  of  the  type  described  and  on  which  Bragg  resonant 
water  waves,  moving  at  the  fluid  velocity  of  the  carrier  wave,  are  superimposed.  For 
simplicity,  it  will  be  assumed  that  there  is  no  interaction  between  the  scatterers  and  the 
carrier  wave,  and  the  scatterers  have  constant  amplitude  throughout.  A  radar  illuminat¬ 
ing  such  a  surface  (Fig.  2)  at  a  depression  angle  y  will  observe  a  doppler  shift  because 
of  the  component  of  fluid  velocity  vK  given  by 


VR 


(D ackz  cos  (V+  y)  +  ua7ke2kz  cos  y 


(6) 


RADAR 


Fig.  2  -  Geometry  of  radar  scattering 
from  a  deterministic  long-crested  com¬ 
posite  surface 


Doppler  spectra  can  be  calculated  by  weighting  the  amplitude  of  a  spectra  of  veloci¬ 
ties,  Eq.  (6),  according  to  the  backscattering  cross  section  given  by  slightly-rough 
theory,  according  to  the  local  angle  of  incidence  of  the  electromagnetic  wave  to  the  water 
wave  13].  Doppler  spectra  calculated  by  this  procedure  have  yielded  broader  spectra  for 
vertical  polarization  than  for  horizontal  polarization,  while  measurements  show  the  con¬ 
verse  Is  true.  The  effect  of  "tilts"  in  the  backscattering  cross  section  (4,5]  would  in¬ 
crease  the  bandwidth  for  the  horizontal  spectra,  and  in  the  limit  it  is  reasonable  to  as¬ 
sume  that  this  deterministic  composite  scattering  model  would  yield  the  same  bandwidth 


*It  is  not  essential  to  use  the  trochoidal  profile,  but  for  a  given  wave  height  and  period  the  trochoid 
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for  both  polarizations.  Thus,  this  model  could  not  account  for  the  polarization  depend¬ 
ence  of  the  spectral  bandwidth  of  radar  sea-echo.  However,  it  does  give  the  bandwidth  of 
the  spectra,  and  the  wave  height  dependence,  and  their  order  of  magnitude.  This  com¬ 
posite  surface  model  will  now  be  carried  over  to  a  surface  like  the  sea. 

The  sea  is  made  up  of  a  continuum  of  water  waves.  Oceanographers  have  found  that 
as  a  first  approximation,  the  sea  can  be  considered  a  Gaussian  random  moving  surface. 
Longuet -Higgins  [14]  derived  many  statistical  properties  for  such  a  surface,  among  them 
the  probability  density  function  for  the  velocity  of  the  zeros  and  specular  points  (points 
on  the  face  of  a  surface  with  a  given  slope).  The  velocity  of  the  zeros  is  equivalent  to 
the  Eulerian  velocity  for  a  single  water  wave,  to  the  first  order  this  velocity  can  be  taken 
as  the  particle  or  fluid  velocity.  Longuet-Higgins’s  probability  density  function  p(v>  for 
the  velocity  of  the  zeros  and  the  specular  points  is 


p(v) 


l_ 

2 


[(v  -  v)2  +  1/  2] 


3/  2 


(7) 


where  v  Is  the  mean  velocity  of  the  fluid  and  v2  relates  to  moments  of  the  water  waves 
spectrum.  For  a  narrow  band  type  of  spectrum,  carrier  water  wave  plus  secondary  wa¬ 
ter  waves  equivalent  to  a  composite  type  scattering  surface,  v 2  in  both  distributions 
simplifies  to 


V 


2 


(8) 


where  n  20  is  the  variance  of  surface  slopes,  noc  the  variance  of  surface  heights,  v  the 
phase  velocity  of  the  carrier  wave,  k  the  wave  number  of  the  carrier  wave,  and  r  the 
group  velocity  (dco/dk).  For  a  gravity  wave,  r  =  v/2.  The  probability  density  function 
for  the  velocity  of  the  zeros  can  be  taken  as  an  approximation  to  the  doppler  ^pectra.  It 
is  not  possible  to  account  explicitly  for  any  scattering  amplitude  weighting  in  this  ap¬ 
proach. 

In  an  alternative  manner  an  average  fluid  velocity  on  the  sea  can  be  calculated  di¬ 
rectly.  Recalling  from  Eq.  (3)  that  the  classical  expression  for  orbital  motion  is 
v,2  =  u,2  +  w,2  =  dj2a2e 2k*  for  a  mean  level  *  =  0,  it  can  be  averaged  in  wave  number' 
space  using  the  energy  spectral  density  of  the  surface,  (P(k).*  Bye  [15]  obtained  drift 
velocity  in  the  sea  using  a  similar  method.  The  averaging  is  performed  in  the  following 
manner:  since  the  variance  of  height  a7  is  related  to  the  area  under  the  spectrum, 


a2  =  f  W(k)  k  dk  ,  (9) 

orbital  motion  co2a2  is  obtained  from 

v,2  =  f  u2  W  (k)  k  dk  ,  (10) 

a  type  of  scalar  average  in  which  phases  and  directions  of  waves  are  not  included,  and 
wherfc  for  gravity  waves  <»2  -  tk. 

In  Eq.  (10)  the  energy  spectrum  for  the  waves  in  the  sea  is  needed,  and  for  conven¬ 
ience  the  spectrum  for  a  fully  developed  sea  will  be  used.  There  are  various  expressions 


*This  method  was  suggested  by  Dr.  J.  W.  Wright  of  the  Naval  Research  Laboratory  in  private 
communication. 
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for  the  '’equilibrium”  range  for  the  energy  spectrum  for  a  fully  developed  sea  [16,17]. 
These  spectra  are  of  the  form 


»T(k)  •**-«,  (11) 

where  a  varies  from  3.75  to  4.5  and  the  constant  B  varies  from  2.5xl0-3  to  7.5xl0’3. 

in  calculating  orbital  velocity  the  lower  and  upper  limits  in  the  integration  will  be 
taken  to  be  gu~2  and  /3  cos  y,  where  g  is  the  acceleration  of  gravity,  U  the  wind  speed, 
p  =  ?.n/\  the  radar  propagation  constant,  and  y  the  radar  depression  angle. 

The  orbital  velocity,  Eq.  (10),  obtained  for  the  different  spectra  is  of  the  form 

v  «c'(iy/7-3/s .  (i2) 

where  the  constant  c '  and  the  range  of  the  exponent  shown  depends  on  the  spectrum  used. 
Equation  (12),  an  important  result,  shows  the  wave  height  dependence  of  orbital  velocity. 
The  width  of  the  radar  spectral  return,  in  units  of  velocity,  should  be  of  th  ^  order  of 
2(v^  )l/2.  Surprisingly  enough,  if  the  fully  developed  sea  assumption  is  used  in  both  the 
deterministic  model  and  Longuet-Higgins  density  function,. we  obtain  a  similar  wave 
height  dependence  as  Eq.  (12),  except  the  value  for  the  constant  depends  on  the  model 
used. 

The  exact  value  of  the  constant  can  be  obtained  if  the  boundary  value  problem  of  the 
electromagnetic  radiation  scattered  by  the  air-sea  interface  is  solved,  or  from  experi¬ 
mental  results. 


EXPERIMENTAL  RESULTS 

Doppler  spectra  have  been  obtained  from  the  phase  data  of  the  Four -Frequency 
Radar  (4FR)  of  the  Naval  Research  Laboratory.  The  complete  system  has  been  described 
in  great  detail  by  Guinard  et  al.  [  18],  thus  only  a  description  of  whatever  is  necessary  to 
clarify  our  analyses  will  be  included  here.  The  4FR  is  a  phase  coherent  pulse  radar 
which  is  capable  of  measuring  the  elements  of  the  scattering  matrix  of  the  radar  return 
at  428  MHz  (P-band),  1228  MHz  (L-band),  4455  MHz  (C-band),  and  8910  MHz  (X-band) 
(Appendix  A).  One  of  the  two  output  channels  of  the  4FR  is  a  pulse -to-pulse  sample  of 
the  cosine  of  the  phase  of  the  clutter  return.  Doppler  spectra  can  be  obtained  from  the 
phase  data  by  the  cosine  Fourier  transform  of  the  correlation  function  of  the  samples,  or 
by  finding  directly  the  power  spectral  estimates  within  a  frequency  interval  by  a  Fourier 
analysis  of  the  record.  The  second  method  was  used  since  it  gave  added  flexibility,  in 
that  degrees  of  freedom,  frequency  resolution,  and  length  of  record  could  be  separated. 
The  resultant  spectrum  calculated  from  the  phase  data  of  the  radar  return  from  an  air¬ 
borne  system  is  the  convolution  of  the  so-called  platform  motion  spectrum  and  the  sea- 
echo  spectrum.  The  platform  motion  spectrum  has  a  mean  doppler  frequency,  fD,  given 
by  the  well  known  expression  for  an  along-track-looking  radar 

?d  -  ~T  cos  y  <  (13) 

A 

where  ve  is  the  ground  speed  and  a  the  radar  wavelength.  In  general,  the  doppler  fre¬ 
quencies  for  the  L,  C,  and  X  bands  are  greater  than  the  radar  pulse  repetition  frequency 
?PRF,  qnd  thus  the  spectra  appear  at  frequency  centers  (aliases)  fciyt  given  by 

f.„,  =  fp  i  n  frpF  .  (14) 
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where  n  is  an  integer.  If  the  spectra  at  the  aliases  appear  '’unfolded”  within  zero  fre¬ 
quency  and  fpRF/ 2  after  detection,  the  spectra  at  the  aliases  are  exact  replicas  of  the 
true  spectra. 

Also,  because  of  the  aircraft  motion,  at  least  two  mechanisms  act  to  broaden  the  re¬ 
sultant  spectra  and  contaminate  the  intrinsic  spectra  of  sea-echo.  These  are  (a)  platform 
motion  (azimuthal  and  elevation  contributions),  and  (b)  decorrelation  because  oi  relative 
motion  and  finite  resolution  cell. 

a.  Platform  motion  —  Because  the  radar  resolution  cell  on  the  sea  is  finite,  and  the 
ground  speed  varies  throughout  the  resolution  cell,  the  resultant  spectrum  is  broadened. 
To  estimate  the  magnitude  of  this  broadening,  Ridenour’s  [19]  expressions  for  an  along- 
track-looking  radar  will  be  used;  they  are: 


and 


l  vsCr 

iv’ 2~r 


sin 


tan 


Ava  =  i  Vtf(A0)J  COS  y 


(15) 


(16) 


where  c  is  the  velocity  of  light  in  vacuo,  r  the  radar  pulse  width,  h  the  aircraft  altitude, 
A0  the  azimuthal  antenna  beamwidth  (half-power),  and  y  the  depression  angle. 

b.  Decorrelation  because  of  relative  motion  and  finite  resolution  cell --This  mecha¬ 
nism  was  discussed  by  Barrick  [20]  and  appears  because  from  one  pulse  to  the  next  the 
resolution  cell  slides  along  the  ground  or  sea  and  cells  for  two  adjacent  pulses  overlap 
only  partially,  introducing  a  decorrelation  of  the  scatterers  responsible  for  the  scatter¬ 
ing.  The  correlation  of  overlapping  scatterers  in  a  pulse -to-pulse  basis  decreases  lin¬ 
early  as  a  function  of  aircraft  velocity  and  time.  The  Fourier  transform  of  such  a  cor¬ 
relation  function  gives  a  spectrum  which  has  an  approximate  half-power  width  Avm 
(velocity  units)  of 

...  Av*  cos  y  fl7) 


Of  course,  atmospheric  turbulence  will  also  affect  the  aircraft  velocity  and  attitude. 
The  antenna  mount  was  a  gyro-controlled  stabilization  system,  and  long-term  attitude 
perturbations  should  be  completely  canceled,  but  without  an  independent  measurement  of 
short  time  perturbations  it  is  not  possible  to  separate  their  magnitude  from  statistical 
fluctuations. 

To  calculate  the  total  broadening  because  of  the  aircraft  motion,  assumptions  of  in¬ 
dependence  between  the  different  mechanisms  and  Gaussian  shaped  spectra  have  been 
made.  Thus,  the  total  broadening  Av  can  be  calculated  from 

Av  ^  [(Av,)3  +  (Av.)»  ♦  (Avm)2]  1/1  .  (18) 

At  grazing  incidence  P-Band  has  the  largest  broadening  and  with  increasing  depression 
angle  all  frequency  bands  attain  the  same  broadening,  Figure  3. 


Bandwidth  of  Radar  Sea-Echo  as  a  Function  of  Sea-State 

Using  the  assumptions  of  independence  and  that  the  shape  of  the  component  spectra 
areGaunsian,  the  nondesirafcle  broadening  can  be  removed  if  its  magnitude  is  small 
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Fig.  3  -  Total  broadening  in  doppler  spectra  because  of 
aircraft  motion  and  finite  resolution  cel'  (aircraft  veloc¬ 
ity:  300  ft/sec,  radar  pulse  width:  0.25  ^sec) 
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compared  to  the  sea-clutter  spectral  width.  To  estimate  the  aircraft's  velocity,  the 
mean  doppler  frequency  from  the  P-band  spectra,  which  is  unambiguous,  has  been  used. 
Any  orbital  motion  shift  will  be  small  compared  to  the  aircraft's  velocity. 


In  Appendix  B  the  sea  conditions  in  which  the  data  were  taken  and  also  other  infor¬ 
mation  pertinent  to  the  spectra  calculations  are  presented.  In  Appendix  C  a  partial  col¬ 
lection  of  the  radar  spectra  used  to  obtain  the  data  has  been  compiled. 


In  Figs.  4  and  5  the  spectra  as  a  function  of  wave  height  for  the  C'-band,  vertical  and 
horizontal  polarizations,  are  shown.  In  the  figures  the  first  letter  in  VV,  HH,  HV,  and 
VH  stands  for  the  polarization  of  the  transmitted  signal,  and  the  second  letter  stands  for 
the  polarization  of  the  received  signal,  where  V  stands  for  vertical  and  H  for  horizontal. 


In  analyzing  the  doppler  spectra,  their  half-power  width  was  measured  directly  from 
the  graphs.  An  actual  calculation  of  the  spectral  variance  to  estimate  width  was  aban¬ 
doned  because  a  direct  measurement  seemed  to  bo  a  fairer  estimate  of  the  actual  spectral 
width  for  the  type  of  data  analyzed.  The  measured  spectral  widths  were  corrected  by 
subtracting  the  square  of  Eq.  (18)  from  the  square  of  the  measured  widths. 


The  broken  line  in  Fig.  6  is  added  to  show  the  wave  height  dependence  H 1  -‘on  orbital 
velocity,  in  agreement  with  Eq.  (12).  The  constant  was  selected  arbitrarily  to  fit  the  data; 
however,  the  scattering  models  predict  the  constant  to  be  anywhere  between  1  and  4. 


In  Fig.  7  the  points  obtained  for  all  the  radar  frequencies  have  been  mciuoeu.  t  .’ 
to  2  ft  waves  the  P-band  data  are  absent  because  platform  motion  was  too  large  to  be  re¬ 
moved. 


Note  the  width  for  the  horizontal  spectra  is  broader  than  the  spectra  for  vertical 
polarization.  Hicks  et  ai.  concluded  that  the  width  for  the  horizontal  spectra  was  greater 
than  that  estimated  from  orbital  motion  alone.  Here  in  agreement  with  Ptdgeon,  we  con¬ 
clude  that  the  spectra  for  vertical  polarization  seem  to  be  attributable  to  orbital  motion 
alone,  while  the  horizontal  spectra  are  broader  and  noisier,  and  should  be  attributed  to 
orbital  motion  and  wind-blown  spray  at  the  air-sea  boundary. 


RELATIVE  AMPLITUDE  (DECIBELS) 
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Fig.  4  -  Doppler  sp€  ra  of  radar  sea-  Fig.  5  -  Doppler  spectra  of  radar  sea- 

echo  broadening  because  of  sea  wave  echo  broadening  because  of  sea  wave 
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Fig=  6  -  Half-power  bandwidth  of  radar  sea- 
echo  spectra  as  a  function  of  significant  wave 
height  of  the  sea  (C-band) 
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Fig.  7  -  Summary  of  corrected  half-pov/er 
bandwidth  of  radar  sea-echo  spectra  from 
the  NRL  4FF  as  a  function  of  significant 
wave  height  of  the  sea  (points  are  averages 
over  depression  angles  between  5  and  30 
degrees) 


HALF  -  POWER  BANDWIDTH  (  FEET  PER  SECOND) 


Another  effect  apparent  in  the  data  is  that  the  width  of  the  spectra  decreases  with 
increasing  radar  frequency;  this  is  surprising,  since  the  converse  seems  to  be  much 
more  reasonable  at  first.  However,  thinking  in  terms  of  the  composite  surface  scatter¬ 
ing  model,  the  effective  rms  slope  of  the  carrier  wave  for  the  P-band  should  be  smaller 
than  that  for  the  higher  frequencies.  Thus,  all  velocities  on  the  surface  would  have  a 
near-constant  amplitude  weight  for  these,  with  the  net  result  of  a  broader  spectra.  For 
calm  sea  condition,  1  to  ft  waves,  the  frequency  dependence  was  weaker  than  in  Fig.  8 
(close  to  unity)  for  all  frequencies  except  the  X-band,  which  had  a  similar  dependence  to 
that  shown  in  the  figure.  The  bandwidth  of  the  spectra  for  vertical  polarization  seems  to 
be  independent  of  depression  angle  for  angles  between  5  and  30  degrees. 

Another  possibility  for  the  broader  spectra  at  the  lower  radar  frequencies  is  the 
longer  length  of  data  required  to  have  a  smaller  frequency  resolution  (Table  B3,  Appen¬ 
dix  B);  any  perturbation  in  the  aircraft  velocity  would  be  much  more  evident  in  the 
spectra  for  these  lower  frequencies.  However,  the  L  and  C-bands  both  use  the  same 
length  of  data  and  the  spec  tra  for  the  lower  frequency  are  still  broader. 


I0Z  I03  10* 

FREQUENCY  (MEGAHERTZ] 


Fig.  8  -  Corrected  half-power  bandwidth  cf 
radar  sea-echo  spectra  as  a  function  of 
radar  frequency  (vertical  polarization) 
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Fig.  9  -  Polarization  dependence  of  the  half¬ 
power  bandwidth  of  radar  sea-echo  spectra 
(21  July  1965) 


In  Fig.  9,  data  from  21  July  1965,  show  the  ratio  of  bandwidth  for  horizontal  polari¬ 
zation  (HH)  to  that  of  vertical  polarization  (W)  decreasing  with  depression  angle.  The 
ratio  of  the  bandwidth  for  the  crosspolarized  spectra  to  that  of  the  vertically  polarized 
spectra  is  smaller  than  the  previous  ratio.  Both  of  these  ratios  seem  to  be  frequency 
dependent.  In  Fig.  10  the  bandwidth  for  horizontally  polarized  spectra  as  a  function  of 
depression  angle  is  shown  explicitly. 


Doppler  Frequency  ~f  Radar  Sea-Echo 

Although  the  bandwidth  of  the  spectra  of  radar  sea-echo  could  not  be  obtained  di¬ 
rectly  and  certain  assumptions  were  necessary  to  recover  it,  the  differential  doppler  be¬ 
tween  the  various  spectra  is  more  exact. 

The  averages  for  the  differential  doppler  as  a  function  of  wind  speed  in  the  C  and  X 
bands  follow  the  general  trend  found  by  Pidgeon  (Fig.  11).  The  important  fact  found  here 
is  that  the  differential  doppler  is  frequency  dependent  as  shown  by  the  figure.  No 
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Fig.  10  -  Half-power  bandwidth  of  radar 
sea-echo  spectra  as  a  function  of  de¬ 
pression  angle  (horizontal  polarization) 
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Fig.  11  -  Average  differential  velocity  of  radar  sea-echo 
spectra  (between  spectra  for  HH  ai;d  VV  polarizations) 
as  a  function  of  wind  speed,  for  depression  angles  smaller 
than  10  degrees 


conclusive  trend  was  found  in  the  data  for  the  days  29  July  1965  and  11  February  1969. 
The  differential  doppler  between  VH  and  W  was  small  for  small  winds  and  increased 
with  wind  speed,  but  its  magnitude  was  still  smaller  than  that  between  (HH  and  W).  As 
expected,  no  noticeable  differential  doppler  was  observed  between  VH  and  HV  for  aU 
wind  speeds  (reciprocity). 


DISCUSSIONS  AND  CONCLUSIONS 

The  experimental  results  on  doppler  spectra  of  radar  sea-echo  from  the  4FR  show 
that  the  mechanisms  responsible  for  sea-clutter  are  quite  complex.  In  particular,  the 
measurements  show  that  several  of  the  parameters  of  the  spectra  are  frequency  and  de¬ 
pression  angle  dependent. 

The  simple  scattering  models  developed  in  terms  of  composite  surfaces  seem  to  be 
most  applicable  to  the  spectra  observed  for  vertical  polarization.  For  example,  the 
measured  bandwidth  of  the  spectra  has  a  wave  height  dependence  similar  to  that  obtained 
from  the  models,  arid  it  is  independent  of  radar  depression  angles  between  5  and  30  de¬ 
grees,  which  is  expected  from  the  models. 

In  the  experimental  results,  the  spectra  for  horizontal  polarization  are  broader  than 
those  for  vertical  polarization,  whire  the  sc&ltering  models  based  on  purely  orbital  mo¬ 
tion  predict  the  converse.  By  allowing  for  "tilt"  effects  in  the  backscattering  cross  sec¬ 
tion  [4,5],  composite  surface  models  at  most  will  yield  the  same  bandwidth  for  the  spectra 
for  both  polarizations.  A  broader  spectra  for  horizontal  polarization  cannot  be  accounted 
for  from  purely  orbital  motion  effects;  thus  there  must  be  an  additional  mechanism  con¬ 
tributing.  Several  characteristics  of  the  spectra  seem  to  indicate  that  spray  over  the 
water  might  be  responsible  for  the  broadening,  e.g.,  the  noisy  appearance  of  the  spectra 
for  horizontal  polarization,  similar  to  the  general  appearance  of  spectra  when  rain  is 
present,  and  the  frequency  dependence  of  the  magnitude  of  the  differential  doppler  be¬ 
tween  HH  and  VV. 
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To  account  for  the  width  of  the  spectra  for  horizontal  polarization,  the  scattering  of 
electromagnetic  waves  from  an  air -spray-sea  interface  should  be  investigated.  As  there 
is  very  little  quantitative  work  published  on  spray  on  the  sea,  this  ?s  an  area  in  which  a 
basic  study  might  be  fruitful. 

In  Ti'ig.  12  the  differential  velocity,  as  a  function  of  depression  angle,  is  shown.  This 
is  the  first  time  that  a  variation  of  differential  velocity  with  depression  angle  has  been 
obtained.  To  arrive  at  a  more  definite  trend,  additional  measurements  in  this  area  would 
be  desirable. 


DEPRESSION  ANGLE  (DEGREES) 

Fig.  12  -  Average  differential  velocity 
of  radar  sea-echo  spectra  (between 
spectra  for  HH  and  VV  polarizations) 
as  a  function  of  depression  angle  (21 
July  1965) 


Figures  13  through  15  illustrate  directly  the  wind  dependence  of  the  differential 
doppler  for  different  radar  frequencies. 

A  summary  of  the  results  obtained  with  the  4FR  phase  data  indicates: 

1.  The  width  of  the  spectra  of  radar  sea-echo  is  frequency  dependent  ,  and  decreases 
with  increasing  radar  frequency. 

2.  The  bandwidth  is  widest  for  horizontal  polarization,  narrowest  for  vertical  polar¬ 
ization,  and  intermediate  for  the  depolarized  spectra. 

3.  The  bandwidth  of  the  spectra  for  all  the  polarizations  increases  with  wave  height. 

4.  The  bandwidth  for  vertical  polarization  is  almost  independent  of  depression  angle 
between  5  and  30  degrees;  the  bandwidth  of  the  spectra  for  horizontal  polarization  de¬ 
creases  toward  the  width  of  the  vertical  spectra  with  increasing  depression  angle. 

5.  Doppler  is  both  polarization  and  frequency  dependent,  greatest  for  the  C  and  X 
bands  and  decreasing  toward  the  P-band. 

6.  The  differential  doppler  is  wind  dependent  for  small  and  intermediate  winds  (0  to 
30  knots)  and  seems  :o  saturate  at  higher  winds  (above  40  knots). 
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Fig.  13  -  Wind  dependence  of  the  differential 
doppler  of  spectra  of  radar  sea-echo  (X-band, 
29  July  1965) 
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Fig.  14  -  Doppler  spectra  of  radar 
sea-echo  for  L-band,  all  polariza¬ 
tions  (29  July  1965) 


Fig.  15  -  Doppler  spectra  of  radar 
sea-echo  for  C-band,  r  11  polariza¬ 
tions  (29  July  1965) 
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7.  The  bandwidth  of  the  spectra  seems  to  be  greatest  upwind,  intermediate  for  cross- 
wind,  and  smallest  for  downwind. 

Work  will  continue  in  the  development  of  a  more  realistic  electromagnetic  scattering 
model  in  which  an  air  plus  wind-blown  spray  and  water  boundary  is  used.  In  order  to 
achieve  this,  it  would  be  very  desirable  to  have  more  quantitative  results  on  spray  over 
the  sea. 
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Appendix  A 

PARAMETERS  OF  THE  FOUR-FREQUENCY  PADAR  SYSTEM 


Table  A1 


Radar 

Fre¬ 

quency 

(MHz) 

Long-Term 
Frequency 
Stability 
(maximum 
parts /day) 

Short-Term 
Frequency 
Stability 
(maximum 
parts/0.01  sec) 

Peak 

Power 

(kW) 

.. 

1 

Average 

Power 

(W) 

Peak 

Power 

Special 

Mode 

(W) 

Inter - 
pulse 
Phase 
Stability 
(maximum 
degrees) 

Output 

Tube 

428 

(P-band) 

1  in  106 

5  in  109  -108 

35-40 

140 

10 

4 

7651 

tetrode 

(RCA) 

1228 

(C-band) 

1  in  106 

5  in  109  -108 

35-40 

140 

10 

4 

7651 

tetrode 

(RCA) 

4455 

(C-band) 

1  in  iO6 

5  in  109-108 

25 

100 

10 

5 

SAC-290 

Klystron 

(Sperry) 

8910 

(X-band) 

1  in  106 

5  in  109  -108 

40 

160 

10 

6 

V24C 

Klystron 

(Varian) 
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Table  A2 

Principal  Parameters  of  Four-Frequency  Radar  System  — 
Signal  and  Display  Parameters 


Repetition  Rates 

Submultiples  of  81.959  kHz: 

100.45,  200.90,  301.35  394.07,  512.30, 

602.71,  683.07,  788.16,  1463,  2926  Hz 

Pulse  Widths 

0.1,  0.25,  0.5,  1,  2  psec 

Noise  Figure 

9-11  including  all  rf  losses 

IF  Amplifier 

Center  Frequency 
Bandwidth 

Type 

Quantity 

37  MHz 

10,  4,  2,  1,  0.5  MHz 

lin-log,  25  and  50  dB  dynamic  range 

Limiting 

Two  —  One  for  horizontally  polarized  signals 
—  One  for  vertically  polarized  signals 

Range  Gate 

Width 

Increment  Readout 

0.1,  0.2,  0.5,  1,  3,  5  psec 

1  yard 

Range  Marks 

Width 

Spacing 

1.0,  2.0,  5.0,  20  ix sec 

1,  5,  20  nautical  miles 

Video  System 

Consoles 

Displays 

CRT  Size 

Total  3:  2  operators  and  1  master  control 
Dual-beam  A  scope,  B  scope  and  PPI  for  each 
console.  All  flat  face  tubes. 

5  inches 

j  Signal  Outputs  (ungated 
and  gated) 

Vertical  polarization  signal  amplitude 

Horizontal  polarization  signal  amplitude 
Vertical/horizontal  polarization  phase 

Horizontal  polarization/transmitted  signal  phase 
Vertical  polarization/transmitted  signal  phase 

Maximum  Video  Amplitude 

2  volts,  peak-positive 

Appendix  B 

SEA  CONDITIONS  AND  PERTINENT  INFORMATION 


Table  B1 

Summary  of  Sea  Conditions 


Measured  RMS 

Visual  Wave 

Wind 

Date 

Wave  Height 

Height 

Speed 

(ft) 

(ft) 

(knots) 

July  29,  1965 

0.39-0.44  (APL) 

1-2 

1-2 

0.40-0.62  (APL) 

2-3 

15 

July  21,  1965 

3-5 

8-11 

'  1 

Dec.  9,  1964 

1.022-1.077  (APL) 

6-8 

12-15 

I  eb.  20,  1969 

16 

29 

Feb.  11,  1969 

26 

46-48 

Table  B2 

Other  Radar  Parameters  During  Measurements 


Run 

Depression 

Angle 

(degrees) 

Altitude 

(ft) 

Antenna 

Azimuth 

(degrees) 

Indicated 

Aircraft 

Speed 

(knots) 

Aircraft 

Heading 

(degrees) 

— 

December  9, 

1964  (PRF  = 

788  Hz,  Pulse  Width  =  0.5  Msec, 

Wind  Direction  =  350  degrees) 

161 

3 

1500 

200  Nominal 

355 

1  162 

5 

1500 

0 

200  Nominal 

170 

!  167 

10 

2000 

0 

200  Nominal 

355 

i  170 

30 

2000 

0 

200  Nominal 

170 

170 

1 

30 

2000 

0 

200  Nomiral 

345 

i 

July  21,  1 

9G5  (PRF  -  788  Hz,  Pulse  Width  0.25  ^sec, 

! 

Wind  Direction  =  100 degrees) 

124A 

10 

5590 

0 

200  Nominal 

270 

124B 

20 

5590 

0 

200  Nomi  d 

270 

124C 

30 

5590 

0 

200  Nominal 

270 

1 24D 

45 

5590 

0 

200  Nominal 

. 

(Table  continue?! 
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Table  B2  (Continued) 


Run 

Depression 
’  ogle 
(degrees) 

Altitude 

(ft) 

Antenna 

Azimuth 

(degrees) 

Indicated 

Aircraft 

Speed 

(knots) 

Aircraft 

Heading 

(degrees) 

125A 

10 

5590 

0 

200  Nominal 

84 

125B 

20 

5590 

0 

200  Nominal 

84 

125C 

30 

5590 

0 

200  Nominal 

84 

125D 

45 

5590 

0 

200  Nominal 

84 

131A 

10 

5590 

0 

200  Nominal 

10 

ldlB 

20 

5590 

0 

200  Nominal 

10 

131C 

30 

5590 

0 

200  Nominal 

10 

131D 

45 

5590 

0 

200  Nominal 

10 

133A 

10 

5590 

0 

200  Nominal 

15 

133B 

20 

5590 

0 

200  Nominal 

15 

133C 

30 

5590 

0 

200  Nominal 

15 

133D 

45 

5590 

0 

200  Nominal 

15 

February  11, 

1959  (tRF  =  683  Hz,  Pulse  Width  =  0.5  jusec 
Wind  Direction  =  290 degrees) 

147 

5 

1500 

0 

195 

308 

150 

10 

1500 

0 

195 

308 

151 

20 

1575 

0 

190 

308 

154 

20 

1575 

180 

190 

308 

155 

10 

1575 

180 

190 

308 

156 

5 

1575 

180 

190 

308 

157 

5 

1575 

0 

190 

128 

160 

10 

1575 

0 

190 

138 

161 

20 

1575 

0 

190 

128 

2 

20 

6750 

350 

192 

220 

3 

30 

6750 

350 

192 

220 

10 

30 

6750 

170 

192 

220 

11 

20 

6750 

170 

192 

220 

36 

45 

12,300 

0 

185 

290 

44 

45 

12,300 

180 

185 

290 

February  20, 

1969  (PRF  =  608  Hz,  Pulse  Width  =  0.5  jusec 
Wind  Direction  =  60  degrees) 

9 

13 

10 

1500 

0 

200  Nominal 

60 

14 

20 

1500 

0 

200  Nominal 

60 

17 

20 

1500 

180 

200  Nominal 

60 

20 

10 

1650 

180 

200  Nominal 

60 

21 

5 

1650 

180 

200  Nominal 

60 

26 

10 

1650 

0 

200  Nominal 

240 

27 

20 

1650 

0 

200  Nominal 

240 

32 

20 

1500 

180 

200  Nominal 

240 

33 

10 

1500 

180 

20C  Nominal 

240 

NOTE: 


0 

0  Antenna  Azimuth 


C  and  X  bands,  forward  looking. 

P  and  L  bands,  backward  looking. 


O 

l',:)  Antenna  Azimuth 


C  and  X  bands,  backward  looking. 
P  and  L  bands,  forward  looking. 
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Appendix  C 

PARTIAL  COLLECTION  OF  SPECTRA 


The  following  graphs  are  in  dB  versus  hertz.  The  labels  near  the  top  of  each  group 
are  of  the  form: 

Date 


Run  number 

Band 

Transmitter 

Receiver 

polarity 

polarity 

Number  of  averages 
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